hep-ph/9907471 
July 1999 



Problems of the Strategy of Regions 



V.A. SmirnovQ 



Nuclear Physics Institute of Moscow State University 
Moscow 119899, Russia 



Abstract 

Problems that arise in the application of general prescriptions of the 
so-called strategy of regions for asymptotic expansions of Feynman in- 
tegrals in various limits of momenta and masses are discussed with the 
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1 Introduction 



When analyzing the leading asymptotic behaviour in various limits a standard strat- 
egy of regions is often used: instead of the integration over the whole space of loop 
momenta only the integration over some specific regions is considered. This strat- 
egy turns out to be possible because the information about the leading behaviour is 
somehow encoded in this integration. A classical example of this procedure is the 
analysis and summation of the leading logarithms in the Sudakov limit [[j]. It was 
argued (and demonstrated for the threshold expansion) in |2| that it is reasonable 
to use the strategy of regions in the most general form and apply it to the whole 
asymptotic expansion, i.e. for any powers and logarithms in an arbitrary limit. In 
this generalized form, the strategy reduces to the following prescriptions: 

(i) Consider various regions of the loop momenta and expand, in every region, the 
integrand in Taylor series with respect to the parameters that are considered small 
in the given region; 

(ii) Integrate the integrand expanded, in every region in its own way, over the 
whole integration domain of the loop momenta; 

(Hi) Put to zero any scaleless integral. 

In the case of off-shell and off-threshold limits, this strategy leads to the well- 
known explicit prescriptions || f|] (see a brief review M) based on the strategy of 
subgraphs. Some subtle points in the application of these general prescriptions of 
strategy of regions to on-shell limits were discussed in ref. |3J. The purpose of this 
paper is to further discuss the status of the strategy of regions using characteristic 
examples of two-loop diagrams. We present two examples of vertex diagrams in 
Section 2. In Section 3, we reformulate the strategy of regions in the language of 
alpha parameters. In conclusion, we discuss the status of the strategy of regions and 
present various recipes for its application. 

2 Two examples 

Our first example is the diagram of Fig. 1 with the masses mi = m 2 = m 3 = m 4 = 
m, m 5 = m 6 = 0, the external momenta p\ = p\ = m2 an d (pi — P2) 2 — ~Q 2 i n the 
limit m/Q — > 0. 

The Feynman integral can be written as 



We use dimensional regularization J7| with d = 4 — 2e. When presenting our results 
we shall omit i7i d ^ 2 per loop and, when writing down separate contributions through 
expansion in e, we shall also omit exp(— 7#e) per loop (je is the Euler constant). 



F 1 (Q,m,e) = J J 



d d kd d l 



(1) 



(I 2 - 2 Pl l)(l 2 - 2p 2 l)(k 2 - 2 Pl k)(k 2 - 2p 2 k)k 2 (k - I) 2 ' 
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Figure 1: (a) Two-loop planar vertex diagram. 



Let us choose, for convenience, the external momenta as follows: 

2 2 
771 7Tb 

Pi = Pi + q^2 , P2 = P2 + qzPi , Pi,2 = (Q/2, =FQ/2, 0, 0) (2) 

so that p 2 = m 2 , pf = 0, 2pip 2 = 2pip 2 = Q 2 - In the given limit, the following 
regions happen to be typical M: 

hard(h): k ~ Q , (3) 

1 - collinear (lc): A; + ~ Q, k_~m 2 /Q, k~m, (4) 

2- collinear (2c): k_ ~ Q, k + ~m 2 /Q, k^m, (5) 

so/if (s): ~ m , (6) 
ultras oft (us): ~ m 2 /Q . (7) 

Here fc± = fco ± fci, k= {k2, k^). We mean by /c ~ Q, etc. that any component of 
is of order Q. 

In the leading order, l/Q 4 , we obtain the following five contributions generated by 
the regions from this list: (h-h), (lc-h), (2c-h), (lc-lc), (2c-2c). (We indicate regions 
for the loop momenta k and I in (|l|) respectively in the first and the second place.) 
In contrast to the diagram with another distribution of the masses and assignments 
of the external momenta, mi = . . . = m 4 = 0, m 5 = m 6 = m ,p\ = p 2 = 0, where 
the last four contributions are not regulated dimensionally so that it is necessary to 
introduce an auxiliary analytic regularization 0, every term can be now considered 
separately, and we symbolically have (2c-h) = (lc-h), (2c-2c) = (lc-lc). 

The (h-h) region generates terms obtained by Taylor expanding the integrand in 
the expansion parameter, m. In the leading order, this is the value of the massless 
planar diagram at p\ = p\ = evaluated in [|1(J : 



c (1) 

c (i) _ c {h-h) 



'(h-h) - (Q2)2+2e 
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Then we have 

d d kd d l 



C 



(i) 

( lc - lci / / ^-2p!/)(/ 2 -2p 2 l)(-2p 1 k)(k 2 -2p 2 k)k 2 (k-l) 2 

d d kd d l 



°(lc-/i) 



(I 2 - 2p 1 l)(l 2 - 2p 2 l)(-2p 1 k)(k 2 - 2p 2 k) 

X k 2 (l 2 -(2p 1 k)(2p 2 l)/Q 2 )- (9) 

Using the technique of alpha parameters and Mellin-Barnes representation we obtain 
the following results: 



C. 



c {1) 

(1) C (lc-lc) 



(lc-lc) - (Q2)2( m 2)2e ' 



a 



c (1) 

(1) _ C (lc-h) 



{lc ~ h) ' {Q 2 ) 2 +"{m 2 Y ' 
c« - -^-^-^_^! + ( e ) (11) 

If we combine all these five contributions we shall obviously obtain a wrong result 
because the poles of the fourth and the third order will not cancel. It turns out that 
the following non-standard region has to be included into the list: 

1-ultracollinear (luc): k + ~ m 2 /Q, fc_ ~ m 4 /Q 3 , k ~ m 3 /Q 2 . (12) 

The missing contributions are (2uc-lc) and (luc-2c) which are equal to each other 
and easily evaluated by alpha parameters for general e. In the leading order, we have 

d d kd d l 



i2 " c J ' ; // (-2pi/)(/ 2 - 2p 2 l){-2p l k){-2p 2 k)k 2 {l 2 - (2p 1 l)(2p 2 k)/Q 2 ) 



c (1) 



(Q 2 ) 2 - 2 ^m 2 )^ ' 



cgL-ic) = -r(e)r(2e)r(3e)r(-4e) 



1 5tt 2 7C(3) 493tt 4 



24 e4 + 48? + ^r + ™r +0(e) - (13) 

Collecting all the seven contributions together we observe that the poles of the 
third and the fourth order in e cancel, and we come to the following result 



(Q 2 ) 2 + 2 ^(Q,m, e ) Q ~°° 
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c m +2c ». (9L)' + 2c m (91) 2 ' + 2c m (91)" 

C (h _ h) + ^C (lc _ h) I 2 I + ^ C (lc-lc) I m 2 I + ^ C (2nc-lc) I 2 I 

, 2 m 2 1 / 5 , , m 2 7r 2 , m 2 >/o .\ 1 

= ln >2?-U ln > + ^ ln ^ +c(3) j^ 

7 , * m 2 An 2 , n m 2 >/ . , m 2 7r 4 ^ . . . . 

+ 8 ln > + -'" V + C(3) 'V + I5 + 0(E) - (14 > 

with the proper coefficient at the double pole which can be evaluated starting from 
the full diagram. 

Note that in the case of the non-planar diagram in the considered limit, with 
Pi = p\ = m2 i there are exactly the same problems with dimensional regularization 
as in the limit with p\ = p\ = ||, |6| . 

The second example is the same diagram Fig. 1 with the masses mi = = 
M, m 2 = m 4 = m, m 5 = m 6 = 0, the external momenta p\ = M 2 , p\ = m 2 
and Q 2 = (pi — p<i) 2 = in the limit m/M — > 0. The corresponding Feynman 
integral F 2 (M, m,e) has literally the same form as ([!]), with the other assignments 
for the external momenta. Let us choose the external momenta as p\ = (M, 0),P2 = 
Mni + ji n 2, with n lj2 = (1/2, =Fl/2, 0, 0). We have 2n 1)2 k = k± for any <i-vector k. 

No 'unusual' regions are relevant here. In the leading order, there are four non- 
zero contributions, corresponding to (h-h), (2c-h),(2c-2c) and (us-2c) regions. (Here 
the same characterization of the regions (3)-(7) in terms of the components k± and 
k, with the substitution Q — > M, is implied.) The (h-h) contribution is obtained by 
expanding the integrand in Taylor series in m 2 . Using again the technique of alpha 
parameters and Mellin-Barnes representation we come to the following result: 



C, 



(2) 

(2) C {h-h) 



~ 126 4 + 12e 2 + 36e + 1440 + ° [ ) ' { &j 



The contribution of the (2c-2c) region is obtained due to the following prescrip- 
tions: 

(a) Expand the propagators 1/(A; 2 — 2p±k) and l/(/ 2 — 2pil) in Taylor series re- 
spectively in k 2 and I 2 , 

(b) Expand each resulting term, which is a function of three kinematical invariants, 
p\ = M 2 , p\ = m 2 and 2pip 2 = m 2 + M 2 , in a Taylor series at the point p\ = and 
2piP2 = M 2 (do not touch p\ 



m 



2\ 



It might seem that we expand in the large mass M but this is just an illusion. 
The (2c-h) contribution is of an intermediate character. In the leading order, we have 

c( 2) r r d d kd d l 



(2c- 2c) 



-2P 2 /)(/ 2 - 2p 2 l)(-2P 2 k)(k 2 - 2p 2 k)k 2 (k - If 
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d&> = 1 1 ^5 m 

(2c ' h) J J (/ 2 - 2 Pl l)(P - 2P 1 l)(-2P 2 k)(k* - 2p 2 k)k 2 (P - (2niZ)2g*) ' 1 ' 

where Pi = Mn\, P 2 = Mn 2 . These contributions are evaluated by the same tech- 
niques as before, with the following results: 

(2) 

c (2) C {2c-2c) 



(M 2 ) 2 (m 2 ) 2e 
_ J 1 ) - 

lc) 
(2) 



J2) _ (1) _ W (17) 

C (2c-2c) — C (2c-2c) — C (lc-lc) ' I 1 'J 



(2) _ C (2c-fe) 



(2c-h) - ( M 2)2+e( m 2) £ ' 

c (2c-h) ~ 8e 4 48e 2 6e 2880 U«J 
The (us-2c) contribution 

c m = f f 

(us - 2c) J J (-2P 2 l)(P-2p 2 l)(-2p 1 k)(-2P 1 k)k 2 {P-(2n 1 k)^) 1 ; 

is evaluated in gamma functions, with a result which is closely related to the (2uc-lc) 
and (luc-2c) contributions in Example 1: 

(2) 

£r(2) _ C (us-2c) 



(us-2c) ( g 2 )2-e( m 2) 3e ' 

(2) . _(1) __(1) 

(us-2c) ~ u (2uc-lc) u (luc-2c) m 



c (2) - c (1) -c (1) (20) 



Collecting all the four contributions together we observe that the poles of the 
third and the fourth order in e cancel, and we come to the following result 

(M 2 ) 2+2e F 2 (M,m,e) M ~°° 

(2) u (2) (M 2 \\ (2) (M 2 \ 2 \ (2) /M*\* 

C (h-h) + C (2c-h) y m 2j + C (2c-2c) y m 2j + C (« S "2c) ^ m 2 J 

, 9 m 2 1 / 1 , , m 2 7r 2 , m 2 _ ,\ 1 
13, 4 m 2 57r 2 l2 m 2 3C(3) , m 2 vr 4 ^. . 

with the proper coefficient at the double pole which can be evaluated starting from 
the full diagram. 
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3 The language of the alpha representation 



There are also contributions of the (h-2c) region (in both examples) and (h-lc) region 
(in the first example) which can be evaluated in gamma functions for general e and 
start from the next-to- leading order, m 2 . Consider now another choice of the loop 
momenta where the momentum I is chosen as the momentum of the 2nd line. Then 
we can recognize a non-zero contribution from the (h-s) region (k is hard and this new 
momentum / is soft). However this is nothing but a double counting because these 
contributions are identical: this can be seen by analyzing this shift of the variables 
/ — > l—p2- This example shows that one has to be careful when testing different choices 
of the loop momenta. Of course, there is no sense of choosing different assignments 
of the loop momenta for the hard regions. But for other types of regions certain 
choices of the loop momenta are crucial — see, e.g., two-loop examples in the case 
of the threshold expansion in ||, in particular, in the case of the ultrasoft regions. 
Here is an another example where a different assignment of the loop momenta is 
important: consider Fig. 1 with mi = . . . = 777,4 = 0, 777.5 — m 6 = Tn,pf = p 2 = 
and m 2 /(pi — P2) 2 — > 0. Then it is necessary to take into account the (h-s) region 
with hard k and the soft momentum of the 6th line || . (It also contributes from the 
next-to- leading order.) 

There is a possibility to completely avoid such double counting by turning to 
the alpha representation and using a similar strategy of regions in the language of 
alpha parameters. This representation for a /i-loop scalar diagram T with powers of 
propagators a\ has the form 



F T (qi 



,Qn) 



:-ir 



ja+h(l-d/2) 7r hd/2 



poo 

x / d«i 
/o 



[°° da L TT a ^ D -d/2 e iA/D-iJ2 mfa t 
JO , 



(22) 



where 



T leT 



A = E 

T 



(Eft) 

l£T \ i / 



(23) 
(24) 



In (p3|), the sum runs over trees of the given graph, i.e. maximal connected subgraphs 
without loops, and, in Q2~4]), over 2-trees, i.e. subgraphs that do not involve loops and 
consist of two connectivity components. The sum of momenta present in (p4|) goes 
over the external momenta that flow into one of the connectivity components of the 
2-tree T. The products of the alpha parameters involved are taken over the lines that 
do not belong to the given tree T. In the non-scalar case, there appear additional 
factors in the integrand of the alpha representation — see, e.g., [IT 
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The strategy of regions in the alpha representation is formulated in the same way 
as in the integrals in the loop momenta. One has to consider alpha parameters to be 
of different order measured in terms of given masses and kinematical invariants. Let 
us illustrate the new language using the two above examples. The function ( p3|) is 
the same in both cases 

D — {a\ + a 2 )(a 3 + a 4 + a 5 ) + a 6 (o;i + a 2 + a 3 + a 4 + a 5 ) . (25) 

The functions (^4j) are 

M = [(«i + ct 3 )(a 2 + a 4 )a 6 + («i + a 2 )a 3 a 4 + aia 2 (a 3 + a 4 + a 5 )]Q 2 

+ [(ai + a 2 )((a 3 + a 4 )(a;i + a 2 + a 3 + a 4 ) 

+ (cti + a 2 )a 5 ) + (ai + a 2 + a 3 + a 4 ) 2 a 6 ]m 2 , 
A 2 = [M 2 («i + a 3 ) + m 2 (a 2 + a 4 )][(«i + a 2 + a 3 + a 4 )a 6 + («i + a 2 )(a 3 + a 4 )] 

+ [M 2 «x + m 2 a 2 ](«i + a 2 )a 5 . (26) 

The contribution of the momentum space (h-h) region is obtained when consid- 
ering all the a\ as being of the same order. (Note that the functions in the alpha 
representation are homogeneous so that it suffices to take into account only the rela- 
tive order of the alpha parameters.) Thus the form D is not expanded and the part 
of the exponent proportional to m 2 is expanded in Taylor series in m. 

In Example 1, we reproduce the contributions of the previously considered mo- 
mentum space regions as follows: 

(lc-h) -> {ct 4i5 ~ m° , ai )2 ,3,6 ~ m2 }\ 

(lc-lc) {« 2 ,4,5,6 ~ rn° , «i,3 ~ m 2 }; 

(h-lc) -> {a 2 ~ m° , ai,3,4,5,6 ~ m 2 }; 

(2uc-lc) — > {a 5 ~ m° , a 3 ~ m 2 , 0:2,4,6 ~ m 4 , a± ~ m 6 }. (27) 

The alpha parameters have dimension mT 2 but to simplify relations we put Q and 
M to one and express all the magnitudes in powers of m. The rest of the regions is 
obtained by permutations. 
In Example 2, we have 

(2c-h) -> {a 4j5 ~ m° , 0(1,2,3,6 ~ rn 2 }; 

(2c-2c) -> {« 2 , 4 ,5,6 ~ "T,° , «i,3 ~ "i 2 }; 

(h-lc) -» {o 2 ~ m° , 01,3,4,5,6 ~ rn 2 }; 

(us-2c) — > {« 5 ~ m° , 02,3,4,6 ~ m 2 , «i ~ m 4 }. (28) 

Note that here we do not have risk to perform double counting of the (h-2c) and 
(h-s) (with another choice of the second loop momentum) regions which was observed 
above. Moreover, the description of the regions in the alpha parametric language is 
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manifestly Lorentz- invariant. However, a visible drawback of this language is rather 
non-trivial description of the (us-2c) and especially (2uc-lc) regions. To find them 
one would need to consider an enormous quantity of various possibilities. 

Let us turn to one more example: the on-shell double box diagram shown in 
Fig. 2. We have pj = 0, i = 1,2,3,4. It was evaluated analytically in ref. [O. An 
explicit analytical algorithm for on-shell double boxes with arbitrary numerators and 
integer powers of propagators was presented in ref. ]n|. In these results were 




3 1 
Figure 2: Double box diagram. 

also checked against first terms of the expansion in t/s, where s = (pi + P2) 2 and 
t = (pi +P3) 2 are Mandelstam variables, with the help of the strategy of regions. The 
external momenta were chosen as 



Pi,2 = (TQ/2, Q/2, 0, 0), r = Pl + p 3 = (-t/Q, 0, yj-t + V/Q\ 0) , (29) 

where s = —Q 2 . The non-zero contributions to the asymptotic expansion in the limit 
t/s — > are (h-h), (lc-lc) and (2c-2c). (The description of the collinear regions is 
given by eqs. (|4]) and @, with the replacement m 2 — > —t.) The (lc-lc) and (2c- 
2c) contributions are not individually regularized by dimensional regularization; a 
natural way to deal with regularized quantities is to introduce and auxiliary analytic 
regularization which is eventually switched off in the sum. 

When analyzing the limit s/t — > 0, it is reasonable to rotate the diagram and 
perform the substitutions t «-> s ,pi — > pi , P2 — > P4, etc. The non-zero contributions 
to the asymptotic expansion in the limit s/t — * are (h-h), (lc-lc), (2c-2c) and (2c- 
lc). Here the introduction of an auxiliary analytic regularization is also reasonable. 
However the corresponding poles in the analytic regularization parameter happen to 
be here up to the second order, as in an example of a non-planar vertex diagram in 
ref. §. 

Let us now describe the contributions to the asymptotic expansion in the limits 
t/s — > and s/t — > in the language of alpha parameters. The functions ( GBp and 
(|24] ) for the double box are 

D3 = (ati + a 2 + a 7 )(a 3 + a 4 + a 5 ) + a 6 («i + a 2 + a 3 + a 4 + a 5 + a 7 ) , (30) 
A 3 = [aia 2 (a 3 + a 4 + a 5 ) + 0:304(0:1 + a 2 + a 7 ) + a 6 (ai + a 3 )(a 2 + a 4 )]s 

+050607^ . (31) 

The (h-h) contributions are reproduced by considering all the alpha parameters of 
the same order. Let us again turn to dimensionless alpha parameters by putting 
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s — — 1 (t — — 1) in the first (second) limit. In the limit t/s —>■ 0, we reproduce 
the contributions of the previously considered collinear momentum space regions as 
follows: 



and, in the limit t/s 



(lc-lc) 


— >■ {"2,4,5,6,7 ~ t° , 


"1,3 ~ 


-*}; 


(2c-2c) 


— * {"1,3,5,6,7 ~ t° , 


"2,4 ~ 


-*} 


-> 0, we 


have 






(lc-lc) 


— >■ {"1,2,3,4,6,7 ~ S° 


, "5 ~ 


-*} 


(2c-2c) 


— >■ {"1,2,3,4,5,6 ~ S° 




-*} 


(2c-lc) 


— ► {"1,2,3,4,5,7 ~ S° 


, a 6 ~ 


-«} 



(32) 



(33) 

We see that, in this example, the description of the collinear contributions in the 
language of alpha parameters is certainly simpler than in the momentum space. 



4 Discussion 

Let us realize that the very word 'region' is understood in the 'physical' sense. In fact, 
it indicates relations between components of the loop momenta expressed in terms 
of the given masses and kinematical invariants. This is clearly not the mathematical 
sense where the region is determined by inequalities. We even do not bother about 
the decomposition of unity, i.e. that our initial integral in the whole space of the loop 
momenta is decomposed into a sum of integrals over all the possible regions which, 
presumably, have zero measure in the intersection of any pair, with their union being 
the whole integration space. 

But the most non-trivial step in the strategy of regions is, probably, the last 
one when all the scaleless integrals are put to zero. Note that this step generally 
does not refer to the use of dimensional regularization. For example, some integrals 
generated by potential contributions within threshold expansion 0] were not at all 
regularized. Still the rule to put such scaleless integrals to zero was experimentally 
checked through examples. 

Up to now there are no counterexamples that would show that the strategy of 
regions generally does not work. Still mathematical proofs in the general case are 
also absent. An indirect confirmation of the strategy is the fact that such proof is 
indeed available |3J] in the case of off-shell and off-threshold limits. (These are limits 
typical for Euclidean space. For this class of regions, it suffices to consider any loop 
momentum to be either hard or soft.) The decisive point of this proof is the analysis 
of convergence of the Feynman integrals. To be more precise, this is the resolution of 
singularities of integrands which is usually performed in the alpha representation. To 
resolve singularities of (|22| ) the whole integration domain is decomposed into sectors. 
These are either a.\ < . . . < a.^ plus sectors obtained by permutations, or more ad- 
vanced sectors associated with one-particle-irreducible subgraphs (and their infrared 
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analogs [H|) of the given graph. After a suitable change of variables in each sector, 
the integrand is factorized and the analysis of convergence reduces to power counting 
in one-dimensional integrals. 

To carry out the analysis of convergence for a given on-shell or threshold limit we 
first need to invent appropriate sectors and sector variables which would provide the 
factorization of the corresponding integrand. At least the sectors mentioned above are 
here insufficient — this can be seen in one-loop examples. It should be also stressed 
that, for each on-shell or threshold limit, this problem of resolution of singularities 
should be solved separately, with sectors and sector variables that are specific for it. 

Suppose that we have to expand a Feynman diagram in some limit typical for 
Minkowski space. The first step is to find all relevant regions that generate non-zero 
contributions. (For the off-shell limits, there is no problem here: we use the well- 
known graph-theoretical language for writing down prescription for expanding the 
given Feynman diagram and just list all relevant subgraphs that are taken from a 
certain family. This task can be even done by computer. In some partial cases, the 
prescriptions for on-shell limits can be also formulated in a graph-theoretical language 
]T5|.) The problem to successfully go through this step is a matter of experience and 



(both physical and mathematical) intuition. The physical flavour of the problem 
is a correspondence of the class of the regions to certain operators and subsequent 
translating the prescriptions for the asymptotic expansion into the operator language. 

When testing various regions it is necessary to be aware of possible double count- 
ing. The combination of this search of the relevant regions both in momentum space 
and alpha integrals looks rather reasonable. But how can we decide that we have 
found all the contributions to the asymptotic expansion? Unfortunately, there is no 
definite answer to this question. At least we can check our results by comparing them 
with one-loop examples where explicit analytical results can be obtained. Sometimes 
comparisons with analytical two-loop results are also available — see, e.g., examples 
in [1^, p. Another important partial check is the cancellation of poles in e, up to 
a certain order, and the analysis of the coefficient at the highest pole which can be 
evaluated by an independent method. 

Thus, if the situation with the poles is unsatisfactory there are at least two options: 

(a) to decide that the strategy of regions breaks down in the given limit; 

(b) to look for missing regions. 

Since the first option has been never realized, it looks reasonable to stay al- 
ways optimistic and continue the search of regions. The regions can be rather 
non-trivial indeed: the ultracollinear region which can be symbolically described as 
luc = (m 2 /Q 2 ) x lc is a strange example. By the way, any 1-ultra-. . . -ultracollinear 
region luc = (m 2 /Q 2 ) h x lc is also a reality for an arbitrary h. Consider the /i-loop 
ladder diagram in the limit of Example 1. Then the contributions of (2u /l_1 c-lu /l_2 c- 
. ..-lc) and (lu /l_1 c-2u ?t ~ 2 c-. . . -2c) regions are non-zero (here h is supposed to be 
even, for definiteness) . 

But imagine now that our terms of expansion satisfy the check of poles. The 
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success is not yet guaranteed because we cannot exclude the existence of a region 
that enters with simple poles in e or even without poles which is insensitive to this 
check. Then one could use numerical checks with numerical evaluation of the initial 
diagram (and, of course, stay optimistic). 

After this advice, let us stress that there is at least an example of the on-shell 
double box diagrams when it turns out easier to evaluate them analytically, rather 



than expand them up to a desired order [[12], |13| . So, the last advice is to try to 



evaluate the diagrams analytically, without any expansion. Still in this case, we can 
use the strategy of regions for crucial checks — see, e.g., [|13 |. 
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